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I.  INTRODUCTION 


In  the  study  of  the  stability  of  a  liquid-filled  projectile,  it  is 
necessary  to  determine  the  flow  field  of  the  liquid  during  the  time  when 
rotational  motion  Is  Imparted  to  it  by  the  container.  This  is  called  the 
spin-up  problem.  Typically,  the  fluid  is  contained  in  a  right-circular 
cylinder.  In  the  present  discussion,  the  fluid  fills  the  cylinder.  The 
coordinate  system  and  some  notation  are  shown  in  Figure  1. 

The  fluid  spins  up  from  rest  until  it  rotates  like  a  solid,  except  for 
the  flow  imparted  by  the  angular  motion  of  the  projectile  which  is  neglected 
in  the  spin-up  problem.  Spin-up  from  rest  is  an  inherently  nonlinear  problem, 
but  the  physics  is  well  understood.  The  basic  work  was  done  by  Wedemeyer.1 
This  and  other  spin-up  problems  are  discussed  by  Greenspan2  and  Benton  and 
Clark.3  More  recently,  the  spin-up  problem  has  been  solved,  without  the 
approximations  of  the  Wedemeyer  model,  using  finite  difference  approximations 
to  the  Navier-Stokes  equations,4 

Wedemeyer' s  model  yielded  a  nonlinear  partial  differential  equation  of 
the  diffusion  type  for  the  azimuthal  velocity,*  V,  as  a  function  of  time  and 
radial  coordinate  but  not  axial  coordinate.  For  large  Reynolds  number,  Re, 
the  diffusion  terms  In  this  equation  can  be  neglected  for  some  purposes  and 
the  solution  for  V  is  elementary;  Wedemeyer  neglected  them  in  his  approximate 
analysis.  For  other  important  applications,  the  diffusion  terms  must  be 
included  giving  a  second-order,  nonlinear,  parabolic  partial  differential 
equation.  The  finite  difference  solutions  of  this  equation,  and  various 
aspects  of  the  theory,  were  presented  by  Sedney  and  Gerber. $ 

Two  of  the  main  appl ications  of  spin-up  theory  are  to  the  study  of  the  free 
oscillations,  the  eigenvalue  problem,  and  forced  oscillations  arising  from  the 


'Definitions  of  symbols  are  given  in  List  of  Symbols  section. 

1.  E.  H.  Wedemeyer,  " The  Unsteady  Flow  Within  a  Spinning  Cylinder ,"  U.S.  Army 
Ballistic  Research  Laboratory ,  Aberdeen  Proving  Ground ,  Maryland,  BRL 
Report  No.  1225,  October  1963  (AD  431846).  Also  Journal  of  Fluid 
Mechanics .  Vol.  20,  Part  3,  1964,  pp.  383-399. 

2.  H.  v.  Greenspan,  The  Theory  of  Rotating  Fluids,  Cambridge  University 
Press,  London  and  New  York,  1968. 

3.  E.  R.  Benton  and  A.  Clark,  Jr.,  "Spin-Up,"  article  in  Annual  Review  of 
Fluid  Mechanics.  Vol.  6,  Annual  Reviews,  Inc.,  Palo  Alto,  California, 

197T. 

4.  C.  W.  Kitchens,  Jr.,  " Navier-Sto«  Solutions  for  Spin-Up  From  Rest  in  a 

Cylindrical  Container,"  U.S.  Army  .allietic  Research  Laboratory,  Aberdeen 
Proving  Ground,  Maryland,  DRL  Technical  Report  ARBRL-TR-02193 ,  September 
1979  (AD  A07711S). 

5.  R.  Sedney  and  N.  Gerber,  " Viscous  Effects  in  the  Wedemeyer  Model  of  Spin- 
Up  From  Rest,"  U .S.  Army  Ballistic  Research  Laboratory ,  Aher  \en  Proving 
Ground,  Maryland,  BRL  Technical  Report  ARBRL  TR-02493 ,  June  1983 

(AD  A12950Q) . 
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angular  motion  of  the  projectile,  the  moment  problem.  The  eigenvalue  problem 
during  spin-up,  i.e,,  determination  of  the  frequencies  and  decay  rates  of  the 
waves  In  the  rotating  f  1  u-* d ,  was  considered  by  Sedney  and  Gerber.6  This 
theory  was  also  used  to  study  the  spin  decay  of  a  liquid-filled  projectile7 
after  ejection  from  the  gun. 

In  these  applications  an  accurate  solution  of  the  spin-up  equation  is 
required.  In  Reference  5  it  was  pointed  out  that  one  of  the  sources  of  error 
in  the  finite  difference  solution  was  the  effect  of  a  discontinuity  in  the 
boundary  data  when  the  rotation  of  the  container  is  imparted  impulsively.  The 
objectives  of  this  report  are  to  provide  the  details  not  included  in  Reference 
5,  extend  the  treatment  of  this  discontinuity,  and  to  give  rules  for  its  prac¬ 
tical  implementation.  The  Wedemeyer  spin-up  model,  the  finite  difference 
method,  and  the  Impulsive  start  assumption  are  discussed  in  detail  in 
Reference  5.  Most  of  the  notation  here  is  the  same  as  in  Reference  5. 


II.  THE  SPIN-UP  EQUATION  AND  BOUNDARY  CONDITIONS 

The  cylinder,  with  radius  a  and  height  2c,  is  filled  with  fluid  and 
initially  at.  rest.  At  time  t  =  0,  it  Is  given  an  angular  velocity  ft  about  its 
axis  which  remains  constant;  i.e.,  the  angular  velocity  history  is  a  step 
function  or  Heaviside  function.  A  description  of  the  resulting  flow  is  re¬ 
quired.  Lengths,  velocities,  pressure,  and  time  are  made  dimensionless  by 
2  2-1 

a,  aft,  pft  a  and  ft  ,  respectively.  In  the  inertial  frame  cylindrical  coor¬ 
dinates  r,  9,  z  are  used  with  velocities  U,  V,  W  in  those  directions,  respect¬ 
ively;  the  origin  of  z  is  at  the  center  of  the  cylinder.  Dimensionless  time 
is  t.  Derivatives  are  Indicated  by  subscripts. 

The  solution  to  the  spin-up  problem  is  governed  by  two  nondimensional 
parameters:  the  aspect  ratio  c/a  and  the  Reynolds  number 


Re  =  ft  a*Vv 

where  v  is  the  kinematic  viscosity  of  the  liquid.  During  spin-up  boundary 
layers,  called  Ekman  layers,  exist  on  the  endwalls.  For  Re  <  105,  approxi¬ 
mately,  the  Ekman  layers  are  laminar,  whereas  they  are  turbulent  for  larger 
Re;  obviously,  this  criterion  is  by  no  means  precise.  The  spin-up  equation 
differs  for  laminar  and  turbulent  Ekman  layers. 


6.  R.  Sedney  and  N.  Gerber,  "Oscillations  of  a  Liquid  in  a  Rotating  Cylin¬ 
der:  Part  II,  Spin-Up,"  U.S.  Army  Ballistic  Research  Laboratory,  Aberdeen 
Proving  Ground ,  Maryland,  BRL  Technical  Report  ARBRL-TR-02489 ,  May  1983 
(AD  A129094). 

7,  C,  W,  Kitchens ,  Jr.,  N.  Gerber,  and  R,  Sedney,  "Spin  Decay  of  Liquid- 
Filled  Projectiles ,"  Journal  of  Spacecraft  and  Rockets ,  Vol.  15,  No.  6, 
Nov-Dec  1978,  pp.  348-314.  See  also  BRL  Report  No.  1996,  July  1977 
(AD  A043275,' 
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It  is  convenient  to  have  a  measure  of  the  time  scale  for  the  spin-up 
process  and  this  depends  on  the  state  of  the  Ekman  layers.  If  the  Ekman 
layers  are  laminar,  the  characteristic  time,  in  seconds,  for  spin-up  is 

fs  =  (2  c/a)  Re1/2/n  (sec) 


or  in  nondlmenslonallzed  form 


ts  =  (2  c/a)  Re1/2. 


If  the  Ekman  layers  are  turbulent,  the  characteristic  time  for  spin-up  is 

£st.  =  (28.6  c/a)  Re1//5/«  (sec) 


or  in  nondimensionalized  form 


t  t  =  (28.6  c/a)  Re1//5. 

These  results  are  derivable  from  lit  aar  spin-up  theory  oi  the  Wedemeyer  solu¬ 
tion  without  diffusion.  The  time  scales  do  not  give  the  time  for  sensible 
conclusion  of  the  spin-up  process;  they  are  a  measure  of  the  time  for  a  change 
of  1/e  from  the  initial  statev 

The  conditions  for  which  the  Wedemeyer  spin-up  model1  is  valid  can  be 
presented  in  terms  of  the  time  scales  involved.  These  are  the  time  for  one 
rotation  of  the  container,  the  spin-up  time  and  the  time  scale  for  diffusion 
across  the  radius  of  the  cylinder.  The  model  was  derived  on  the  basis  of  an 
Impulsively  started  cylindrical  container,  a  condition  which  can  only  be 
approximated  in  practice.  Additional  time  scales  must  be  introduced  to  assess 
the  approximation  to  an  in^  lslve  start;  these  are  discussed  in  Referei .ce  5. 
The  acceleration  of  a  projectile  in  a  gun  tube  is  a  relatively  good  approxi¬ 
mation  to  an  impulsive  start.  Of  course,  the  actual  time  history  of  projec¬ 
tile  angular  velocity  in  the  gun  tube  could  be  used  in  the  spin-up  problem, 
after  suitable  modification  of  the  Wedemeyer  imp  Isive  start  model,  but  the 
finite  difference  calculation  would  require  time  steps  which  would  be  small 
compared  to  the  acceleration  timt,  which  is  typically  0.02  sec,  and  would 
require  a  new  calculation  for  each  gun  tube  twist  and  for  each  zone  charge. 
Such  an  approach  is  less  practical  and  more  time  consummirig  than  using  the 
impulsive  start  approximation. 

For  an  impulsive  start,  Wedemeyer ' s  model1  determines  the  "core  flow," 
i.e.,  the  flow  exterior  to  the  Ekman  layers  and  a  sidewall  (Stewartson) 
layer.  Although  the  action  of  the  endwalls,  through  the  Ekman  layers,  is 
essential  to  the  spin-up  process,  only  the  solution  to  the  core  flow  is 
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determined  by  the  model.  Usi.:g  order-of-magnitude  arguments,  Wedemeyer  showed 
that  the  Navier-Stokes  equations  reduce,  for  the  core  flow,  to 

Vt  +  U  (vr  ♦  V/r)  =  Re-1  [Vrr  +  (V/p) r3  (2.1) 

and 

uz  '  vz  =  pz  '  0  (2.2) 

plus 

Pp  =  V2/r  and  (rU)r  +  r  Wz  =  0. 

For  Re  *  *  he  proposed  neglecting  the  diffusion  terms  in  (2.1)  so  that 


\  +  Uw  <Vwr  *  V>  ■  0  (2-3) 

where  the  subscript  w  denotes  this  approximation. 

To  solve  (2.1)  or  (2.3),  a  relationship  between  U  and  V  is  necessary. 

This  is  called  the  Ekman  compatibility  condition  because  the  Ekman  layer 
suction  must  be  made  compatible  with  the  core  flow.  A  phenomenological 
approach  was  necessary  at  this  point  in  the  theory.  For  laminar  Ekman  layers, 
Wedemeyer  proposed 


U  =  -k4  (r  -  V) 

k£  -  k  (a/c)  Re"1/4  =  2</ts 


(2.4) 


with  k  =  0.443;  k  =  0.5  often  gives  results  in  better  agreement  with  numerical 
solutions  to  the  Navier-Stokes  and  will  usually  be  used  here. 


For  turbulent  Ekman  layers 

U  =  -kt  (r  -  V)8/5 

kt  =  0.035  (a/c)  Re"1/5  *  l/tst. 


(2.5) 


Using  (2.4),  (2.3)  can  be  solved  for  V  (r,  t)  with  V  (r,  0) 
( 1 ,  t)  -  1 ; 


and  V 


2k  {t  2k ,t 

Vw  -  (re  A  -  l/r)/(e  1  -  1) 


for  r  >  e 


-k„t 


=  0 


for  r  <  e 


-k„t 


(2.6) 
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Thus  r  -  pf  = 


separates  rotating  and  nonrotating  fluid  and  Is  called 


the  front.  Using  (2.5)  In  (2.3),  a  numerical  Integration  is  necessary  but  the 
character  of  the  solution  is  the  same  as  when  (2.4)  is  used.  The  solution 
(pp.  8-16  of  Reference  8)  Is 


ktt  =  f1  "  .V/0"8/5  (2.7) 

where  r 

y  =  r  Vw  (y  kept  constant  for  the  integration). 

The  solution  is  completed  by  obtaining  U  from  (2.4)  or  (2.5)  and  W  from  the 
continuity  equation.  The  solutions  (2.6)  and  (2.7)  have  a  shear  discontinu¬ 
ity,  i.e.,  discontinuous  Vw^,  which  would  be  smoothed  out  if  the  diffusion 

terms  were  retained  in  (2.1). 

It  follows  easily  from  (2.1)  and  (2.4)  that 


V  =  f  (r ,  k jt,  Re); 

therefore  k^t  (=  t/t$  for  k  =  0.5)  is  the  natural  time  scale.  For  the 
turbulent  case  k ^  is  replaced  by  kt. 


In 

V(r,  t) 
problem 

accordance  with  (2.2),  V  is  independent  of  z.  To 
initial  and  boundary  conditions  must  be  applied, 
these  are: 

solve  (2.1)  for 

For  the  present 

V  =  0 

for  t  <  0,  0  <  r  <  1 

(2.8) 

V  -  1 

for  r  =  1  t  >  0 

(2.9) 

V  =  0 

f  jr  r  -  0  t  >  0. 

(2.10) 

The  analogy  between  (2.1)  and  the  heat  conduction  equation  allows  on’i  to 
deduce  that  these  are  necessary  and  sufficient  for  determining  a  solution. 
The  form  of  (2.9)  follows  from  the  Impulsive  start  condition.  At  r  =  1,  V 
would  be  a  specified  function  of  t  for  a  nonlmpulslve  start.  Conditions 
(2.8)  and  (2.9)  require  a  discontinuity  in  V  at  the  point  r  =  1,  t  =  0. 


8,  Engineering  Design  Handbook ,  Liquid-Filled  Projectile  Design ,  AMC  Pamphlet 
706-165 3  April  1969  (AD  85S?TUy. 
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A  finite  difference  scheme  for  solving  (2.1)  was  presented  in  Referent 
5.  The  value  of  V  (1,0)  Is  required  but  since  0  <  V  (1,0)  <  1  the  value  is 

undetermined.  If,  for  example,  a  value  of  V  (1,0)  =  0,  1/2,  or  1  is  used. 

there  will  be  a  significant  error  In  V;  the  error  Is  largest  near  r  -  1. 
Because  of  the  diffusive  nature  of  (2.1),  this  error  will  decrease  as  t  in¬ 
creases.  However,  it  was  found  that  the  error  was  often  significant  even 
for  t»  t$.  Obtaining  a  finite  difference  solution  when  a  discontinuity  In 

boundary  conditions  exists  can  always  be  expected  to  require  special  treat¬ 
ment.  The  approach  adopted  here  and  presented  in  Section  IV  is  to  obtain  a 
local  solution  which  includes  the  effects  of  the  discontinuity;  that  solution 
Is  used  to  generate  a  new  Initial  condition.  3efore  doing  that  we  digress  to 
consider  the  manner  In  which  the  discontinuity  Is  treated  for  the  non- 
dlffuslve  equation. 


III.  THE  NONDIFFUSIVE  EQUATION 

The  approximation  of  (2.1)  by  (2.3)  Is  clearly  a  singular  perturbation 
since  the  second  order  equation  Is  reduced  to  a  first  order  equation.  The 
boundary  conditions  (2.8)  and  (2.9)  are  applied  to  (2.3);  (2.10)  cannot  be 
Imposed  but  the  solution  (2.6)  satisfies  it.  The  discontinuity  can  be 
accommodated  in  solving  (2.3),  as  will  be  shown. 

It  is  convenient  to  Introduce  the  circulation,  r  =  r  V  ,  and  the  time 
scale  w 


f  -  t. 

Only  the  laminar  case  will  be  considered  here.  Then  (2.3)  becomes,  using 
(2.4), 


rt.  ••  (r  -  r/r)  rr  =  o . 


(3.1) 


This  first  order  equation  has  a  one-parameter  family  of  characteristics  and 
can  be  solved  bv  the  standard  method  of  characteristics.  These  characteris¬ 
tics  contain  some  information  related  to  the  properties  of  the  solution  of  tne 
original,  diffusive  equation,  (2.1),  and  they  are  called  the  sub-characteris¬ 
tics  of  t.ie  original  equation.9  Along  a  characteristic  of  (3.1)  r  =  r  = 
constant,  0  <  r  <  1,  and  the  equation  of  the  characteristics  is  J 

r2  =  r0  +  (y2  -  rQ)  e"2t’.  (3.2) 

where  y  is  a  constant,  0  <  y  <  1.  This  equation  can  be  written  in  two  parts: 


9. 


J,  D.  Cole ,  Perturbation  Methods  in  Applied  Mathematics, 
Publishing  Company ,  Waltham ,  MA3  1968 . 


Blaisdell 
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(3.3) 


i-  =  f  e_t' .  r0  -  o 

which  curves  issue  from  t'  =  0,  0  <  r  <  1  and 

r2  =  rQ  +  (1  -  rQ)  e“2t',  y  =  1  (3.4) 

which  curves  emanate  from  the  point  r  =  1,  t‘  =0  where  the  discontinuity 
exists.  The  characteristic  common  to  (3.3)  and  (3.4)  is  the  front 


Some  examples  of  the  characteristics  are  shown  in  Figure  2.  For  the  charac¬ 
teristics  to  the  left  of  the  front,  (3.3),  *  -*•  0  as  t'  -*•  »  whereas,  for  those 

to  the  right,  (3.4),  r  I^1/2  as  t‘  «.  Some  discussion  of  the  method  of 

characteristics  was  given  by  Wedemeyer8  and  Weidman.10 

The  set  of  characteristics,  (3.4),  that  originates  at  the  point,  r  =  1, 
t'  =  0  is  called  a  fan.  A  more  familiar  example  of  a  fan  occurs  in  supersonic 
flow  over  a  corner,  the  Prandtl -Meyer  expansion.  All  the  given  data, 

0  <  rQ  <  1,  at  the  discontinuity  is  convected  into  t*  >  0  along  the  fan.  The 

initial  condition  rQ  =  0  is  convected  along  (3.3). 

Using  the  method  of  characteristics  provides  some  insight  into  the 
solution  of  (2.3)  and,  in  particular,  how  the  discontinuity  is  accommodated. 

Of  course,  all  the  information  is  contained  in  (2.6).  It  could  be  used  to 
obtain  an  initial  condition  for  (2.1)  to  treat  the  discontinuity.  (See  Section 
V.)  It  would  be  better  than  assuming  a  value  for  V( 1 ,0)  but  considerably  less 
accurate  than  the  method  discussed  next. 


IV.  THE  SOLUTION  FOR  SMALL  TIME 


Introducing  (2.4)  and  the  time  scale  t'  =  k^t  into  (2.1)  gives 

Vt.  +  (V  -  r)  (Vr  +  V/r)  =  e  [Vrr  +  (V/r)p]  (4.1) 


10.  P.  D.  Weidman,  "On  the  Spin-Up  and  Spin-Down  of  a  Rotating  Fluid:  Part 
1.  Expending  the  Wedeme Model,"  Journal  of  Fluid  Mechanics,  Vol.  77 , 
Part  4,  1976,  pp.  685-7ud. 
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where  e=  1/k^Re.  The  boundary  and  ir.itial  conditions  are  still  (2.8)- 

(2.10).  The  laminar  Ekman  layer  case  will  be  treated  in  detail;  the  results 
for  the  turbulent  case  are  merely  indicated. 

The  domain  of  integration  for  (4.1)  is  shown  in  Figure  3a  together  with 
the  boundary  and  initial  conditions;  the  discontinuity  is  at  point  P: 
r  =  1,  t  =  0.  A  solution  is  sought  which  satisfies  (4.1),  (2.8),  and  (2.9)  in 
the  neighborhood  of  P.  The  approach  is  to  stretch  the  neighborhood  of  P,  by 
means  of  a  singular  transformation,  so  that  the  behavior  of  V  can  be 
examined.  The  coordi nates  (r,  t')  are  transformed  to  (R,  T): 


R  -  (1  -  r)/2  (e  t')1/2 
T  =  t' 


(4.2) 


This  definition  of  R  is  1/2  e^2  that  of  the  R  in  Reference  5.  The  domain  of 
integration  in  the  transformed  plane  is  shown  in  Figure  3b  with  the  trans¬ 
formed  boundaries  of  the  (r,  t')  region  indicated.  The  point  P  transforms 
into  the  positive  R  axis,  the  initial  line  t *  =  0  transforms  to  the  point  at 
infinity,  and  r  =  0  transforms  into  a  hyperbola-like  curve. 


Equation  (4.1)  is  transformed  to  (R,  T)  and  a  solution  near  T  =  0  is 
sought.  The  form  of  the  equation  suggests  the  expansion 


V  -  VQ  (R)  +  T1/2  Vj  (R)  +  ...  (4.3) 

where  VQ  and  Vj  satisfy 

Vq  +  2  R  Vq  =  0  (4.4) 

V'|  +  2  R  V|  -  2  Vx  =  2  V;  [e1/2  +  e"1/2  (1  -  VQ)]  (4.5) 


with  '  =  d/dR.  The  series  (4.3)  is  a  singular  perturbation  as  evidenced  by 
the  fact  that  the  partial  differential  equation  has  been  replaced  by  ordinary 
differential  equations  and,  as  will  be  shown,  the  boundary  condition  at  r  =  0 
cannot  be  satisfied.  In  the  terminology  of  Van  Dyke11  (4.3)  is  a  direct 
coordinate  expansion.  Such  expansions  are  usually  nonuniform,  e.g.,  diverging 
for  large  T  which  is  true  for  (4.3).  For  T  small  enough  the  two-term  series 
(4.3)  gives  an  adequate  approximation  to  V. 

The  boundary  conditions  for  V0  are  obtained  from  (2.8)  and  (2.9): 

V0  (0)  =  1  VQ  (»)  =  0. 


11.  M.  Van  Dyke.  Perturbation  Methods  in  Fluid  Mechanics ,  Academic  Press ,  New 
York ,  NY,  19641 
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The  solution  to  (4.4)  is 


VQ  =  erfc  R. 


(4.6) 


For  r  «  1,  this  is  the  same  as  the  solution  for  the  impulsively  started, 
infinite  plate  (the  Rayleigh  problem)  which  is  to  be  expected  since  curvature 
effects  are  negligible  for  r  »  1.  The  boundary  conditions  for  are 


V*  (0)  =  Vj_  H  =  0 


and  the  solution  to  (4.5)  is  then 


/j  =  e1/2  R  (1  -  erf  R)  +  e“1//2  [R  (1  -  erf2R)  -  iT1/2  e_!r  erf  R].  (4.7) 


The  one-term  solution,  V  =  V0  =  erfc  R  has  the  property  that  V  =  constant 

along  the  parabolas  in  the  (r,  t)  plane,  R  =  constant.  This  is  somewhat 
analogous  to  the  property  of  the  nondiffusive  equation  that  r  =  constant 
along  the  fan  (3.4). 

The  original  objective  was  to  obtain  a  solution  to  (4.1)  in  the 
neighborhood  of  P  but  (4.3),  (4.6),  and  (4.7)  predict  V  for  all  r.  For 

r  =  0,  or  R  *  1/2  (e  T)1/2,  it  is  clear  that  neither  VQ  nor  Vj  is  zero  so  that 

(2.10)  is  not  satisfied.  In  principle  a  solution  near  r  =  0  should  be 
obtained  and  matched  to  (4.3).  In  practice,  for  application  to  the  finite 
difference  problem  this  is  not  necessary,  as  shown  in  the  next  section. 

Plots  of  Vj  for  three  values  of  e  and  VQ  are  presented  in  Figure  4.  Both 
VQ  and  decrease  exponentially  as  R  +  ».  Their  asymptotic  forms  are 

R2 

V0  ~ ( l//i)  e  K  /R  (4.8) 


2 

Vj  ~(l//i)  (e1/2  +  e_1/2)  e"R  .  (4.9) 


For  estimation  purposes  (4.8)  is  in  error  by  1  in  the  second  decimal  place  for 
R  =  1.34  and  by  1  in  the  sixth  decimal  place  for  R  =  3.12;  for  a  moderately 

small  e  =  .0313,  corresponding  to  Re  =  4  x  104  and  c/a  =  3.120,  (4.9)  is  in 
error  by  1  in  the  first  decimal  place  for  R  =  1.34  and  by  11  in  the  fifth 
decimal  place  for  R  =  3.12. 
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V.  APPLICATION  OF  THE  ELUTION  FOR  SMALL  TIME 


If  the  Initial  condition  (2.8)  Is  applied  at  t  =  0.  the  difficulties 
discussed  In  Section  II  arise.  The  limit  of  (4.3)  as  t  or  t  +  0  is  (2.8). 
Therefore,  we  must  apply  (4.3)  at  a  small  time,  the  natural  choice  being  on. 
time  step.  At,  in  the  finite  difference  method.  Although  we  shall  make  that 
choice,  the  method  Is  not  restricted  to  it. 

Either  the  one-term  or  two-term  solutions  in  (4.3)  can  be  used.  We  set 


v  =  v0  (R) 

(5.1) 

V  =  V0  (R)  +  T1/2  Vj  (R) 

(5.2) 

at  t  =  At  or  T  =  At.  Expressing  R  In  terms  oi  the  original  variables  gives 

R  -  (1  -  r)  (Re/At)1/2/2. 


The  functions  of  r  computed  from  either  (5.1)  or  (5.2)  provide  the  initial 
conditions  at  t  =  At.  The  fact  that  V  *  0  at  r  =  0  is  of  no  practical 
concern.  Only  an  approximate  Initial  condition  can  be  expected  and  the 

approximation  at  r  *  0  is  quite  good.  For  example,  if  R  >  3.5,  V  <  10 *6  if 
(5.1)  is  used;  for  a  range  of  At  and  e  the  same  is  true  if  (5.2)  is  used.  For 
r  =  0  this  condition  on  R  is 


(Re/At)1/2  >  7. 


(5.3) 


For  At  =  1,  the  lower  bound  on  Re  is  49.  For  smaller  Re,  At  must  be 
decreased. 

There  is  an  indepenr'  t  test  of  the  permissible  At  in  the  iterative 
finite  difference  calculation.  The  At  is  required  to  be  small  enough  so  that 
no  more  than  4  iterations  are  required  for  convergence  at  any  time  step. 
Experience  has  shown  that  if  the  £  from  the  criterion  is  satisfied,  (5.3)  is 
also.  The  requirement  on  Ar  for  the  finite  difference  calculation  is  based 
mainly  on  experience  and/or  tests;  see  Reference  5  for  a  discussion  of  this. 
There  is  obviously  an  interaction  between  the  Ar  error  in  the  finite  diff¬ 
erence  scheme  and  the  accuracy  of  the  initial  condition.  For  small  enough 
At,  (5.2)  is  more  accurate  than  (5.1),  and  using  it  gives  a  more  accurate  V 
for  a  given  at;  alternatively,  for  a  given  accuracy,  a  larger  Ar 
can  be  used  with  (5.2). 

The  accuracy  of  the  two-term  solution,  for  small  t,  can  be  illustrated  by 
comparing  it  with  the  finite  difference  solution  of  (2.1)  using  (2.4);  a 
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criterion  for  small  t  Is  given  In  the  next  paragraph.  For  Case  1,  with  Re  = 

4  x  10^,  c/a  =  3.120,  and  e  =  .0313,  the  solution  V  (r,  t)  Is  shown  In  Fig,  - 

5  for  .84  <  r  <  1.0  at  t  »  40  with  initial  condition  (5.2)  at  t  =  At  =  5. 

The  solutions  VQ  and  VQ  +  T*^2  Vj  at  t  =  40  are  plotted  and  show  that  the  two- 

term  approximation  Is  still  accurate  at  t  =  40,  The  nondiffuslve  solution, 
(2.6),  Is  also  plotted  to  lllustate  the  remark  at  the  end  of  Section  III. 

In  the  discussion  after  (4.3)  the  likelihood  that  it  would  diverge  for 
increasing  T  was  mentioned.  This  Is  illustrated  in  Figure  6.  The  finite 
difference  result  for  V,  using  (5.2)  at  At  s  5,  Is  shown  vs  t  for  Case  1, 

c 

r  =  .94  and  Case  2,  r  *  .995.  Case  2  parameters  are  Re  =  4  x  10  ,  c/a  =  1, 

and  e  =  .00316.  The  one-term  and  two-term  solutions  are  also  shown.  For  Case 

1/2 

1,  Vn  +  T  V,  is  a  good  approximation  to  V  but  is  beginning  to  diverge  from 

it.  For  Case  2,  the  divergence  is  quite  evident  and  in  fact  VQ  +  T  '  >  1 

for  t  >  42.  Since  V  >  1  is  impossible,  the  two-term  solution  is  not  valid 
near  t"=  42.  An  approximate  expression  for  the  time  of  breakdown  can  be 
obtained  from  the  condition  Vr  (r  =  1)  =  0.  Using  the  two-term  approximation 

for  V  this  gives 


+  T 


1/2 


vl)rlr 


=  C. 


Using  (4.6)  and  (4.7)  this  yields 


Ti/2  [el/2  +  e-l/2 


(1  -  2/tt)]  =  2iT1/2 


or,  for  e  small , 


T  *  10e 


(5.4) 


as  the  approximate  time  of  breakdown.  For  Case  1  in  Figure  6,  this  gives 
t  »  400  and  for  Case  2  t  *  40.  If  the  solution  is  val id  up  to  a  time  which  is 
some  fraction,  C,  of  T  in  (5.4) 

t  *  (10  C/k2)  (c/a)2 


as  a  practical  upper-limit  for  the  validity  of  the  two-term  approximation. 

For  the  cases  considered  C  =  0.5  is  suitable. 

The  crucial  tests  of  the  application  of  (5.1)  or  (5.2)  is  the  accuracy  of 
the  finite  difference  solution  for  times  comparable  to  t$.  Results  will  be 
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shown  for  Case  3  with  Re  =  4974,  c/a  =  3.30,  e  =  .1056,  and  x  =  .443  giving 

=  525.  The  comparison  will  be  given  first  in  terms  of  the  vorticity  c  = 

(r  V)r/2r.  This  is  an  Important  quantity  because  it  occurs,  along  with  V.  .. 

the  perturbation  equations.  Tor  comparison  purposes,  it  has  the  advantage 
that  an  exact  result  Is  known:  cr  (r  =  1)  =  0  fcr  all  t.  This  follows 

directly  from  (4.1).  In  Figure  7  5  vs  r  Is  shown  for  Case  3  at  t  -  100  and 
600;  c  is  obtained  from  the  finite  difference  solution  using  (5.1)  as  initial 
condition  at  at  *  5  and  using  V  (1,0)  =  1.  For  t  =  100  the  errors  resulting 
from  use  of  V  (1,0)  «  1  are  quite  large;  for  t  =  600  the  errors  are  smaller 
but  would  be  unacceptable  In  the  solution  of  the  eigenvalue  problem.  For  t  = 
100,  (r  =  1)  =  -240  for  the  V  (1,0)  =  1  curve;  c(,  (r  =  1)  =  -1.5  for  the 

curve  using  (5.1).  The  results  using  (5.2)  are  not  shown  in  Figure  7 
but  cr  (r  =  1)  =  0.32  is  obtained. 

For  t  =  100,  the  assumption  V  (1,0)  =  1  causes  errors  in  V  of  the  order 
of  10%;  whereas,  for  t  =  600  the  errors  are  about  3%.  Experience  has  shown 
that  errors  of  3  -  4%  in  V  give  erroneous  results  for  the  eigenvalue. 


VI.  THE  TURBULENT  CASE 

The  same  kind  of  approach  can  be  used  to  derive  the  initial  condition  for 
the  turbulent  case.  In  fact  the  one-term  solution  is  the  same  as  for  the 
laminar  case,  (5.1).  Only  the  essential  steps  will  be  given  here.  For  the 
turbulent  case,  (2.1)  with  (2.5)  must  be  solved  for  V  with  the  same  boundary 
conditions  (2.8)  -  (2.10).  Analogous  to  (4.1),  V  satisfies 


V  -  (r  -  V)8/5  (Vr  +  V/r)  *  et  [Vpr  +  (V/r)p]  (6.1) 

where  now  t'  =  ktt  and  =  l/kt  Re.  Again  the  transformation  (4.2)  and  the 
expansion  (4.3)  are  used.  The  equation  for  VQ  and  its  solution  are  the  same 
as  (4.4)  and  (4.6),  respectively.  The  equation  for  Vj  is  now 

VJ  +  2  R  V'  -  2  =  2  V;  [ej/2  +  e'l/2  (1  -  VQ)8/5]  (6.2) 


With  V0  =  erfc  R  It  is  unlikely  that  a  closed  form  solution  to  (6.2)  can  be 

obtained  as  was  done  for  the  laminar  case;  the  homogeneous  solution  of  (6.2) 
is  the  same  as  before.  A  numerical  integration  of  (6.2)  is  certainly  feasible 
but  has  not  been  incorporated  into  the  Initial  condition.  For  the  turbulent 
case  only  the  one-tern,-  solution,  V0,  has  been  used. 


VII.  CONCLUSIONS 


For  an  impulsive  start,  the  boundary  conditions  for  the  spin-up  equation 
prescribe  a  discontinuity  in  the  velocity  on  the  boundary.  If  this  discon¬ 
tinuity  is  not  treated  properly,  the  finite  difference  solution  of  the  spin-up 
equation  will  be  in  error.  Although  the  errors  continually  decrease  (because 
the  spin-up  equation  is  of  the  diffusion  type)  they  may  be  significant  e^en 
for  one  spin-up  time.  The  treatment  of  the  discontinuity  proposed  here  con¬ 
sists  of  determining  a  local  solution  near  the  discontinuity  and  using  that  as 
an  approximate  initial  condition.  The  degree  of  approximation  has  been  in¬ 
vestigated  and  limits  of  applicability  determined.  Tests  of  this  approximate 
initial  condition,  over  a  range  of  Re  and  c/a,  have  shown  that  the  errors  in 
the  solution  are  negligible  at  about  one-half  of  the  spin-up  time;  the  time 
for  negligible  errors  can  be  decreased  by  adjusting  the  parameters  of  the 
initial  solution. 


19 


-  Finite  Diff.,  using  (5.?) 

.  V0+(T1/2)  Vj 

-  Nondlffuslve  V 


Figure  5.  Comparison  of  the  Finite  Difference  Solution  for  V  at 
t  =  40  With  the  One-Term  and  Two-Term  Solutions,  (5.1) 
and  (5.2),  and  the  Nondiffusive  Solution  (2.6). 

Re  s  4  x  104,  c/a  =  3.120,  e  =  0.0313. 


Case  2 

r=.995 

Finite  Difference  V  with 

(5.2)  at  At  =  5.0 

Case  1 

. VTl/2  V1 

r=.940 

1  f 

- v0 

Figure  6.  Finite  Difference  Solution  V  vs  t  for  Cases  1  and  2  Com¬ 
pared  With  the  One-Term  and  Two-Term  Approximate  Solutions. 

Case  1:  Re  =  4  x  104,  c/a  =  3.120; 

Case  2;  Re  =  4  x  105,  c/a  =  1. 
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LIST  OF  SYMBOLS 


cross-sectional  radius  of  cylinder 
half-height  of  cylinder 

s  k  (a/c)  Re'1/2  (See  (2.4).) 

i  0.035  (a/c)  Re"1/5  (See  (2.5).) 

number  of  subintervals  in  r  in  finite-differt  ice  solution 


pressure/( pa^fr)  in  Wedemeyer  model  (See  (2.2),) 
radial  coordinate/a 

location  jf  front  in  Wedemeyer  nondiffusive  ir'del 


i  (1  -  r)/(4et* )1/2 


=  a^n/v,  Reynolds  number 
time  x  Q 


fst  n 


=  k£t  (or  k^t) 


*  * 


J,  v,  w 


uw,  vw 


=  (2  c/a)  Re^2/Q,  characteristic  spin-up  time  for  laminar 
Ekman  layer 

=  (28.6  c/a)  Re1/5/fl,  cnaracteristic  spin-up  time  for 
turbulent  Ekman  layer 

=  t1  (see  (4.2)) 

radijl,  azimuthal,  axial  velocity  components  *  l/(aft)  of 
Wedemeyer  model  spin-up  flow  with  diffusion  (See  (2.1)  and 
(2.2).) 

radial  and  azimuthal  velocity  components  x  l/(an)  of 
Wedemeyer  model  spin-up  flow  without  diffusion  (See  (2.3). 
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V0(R),  V^(R)  zeroth  and  first  order  coefficients  of  T1^2  in  (4.3) 

z  axial  coordinate  (z  =  0  at  cylinder  midplane.) 

y  constant  in  characteristic  equation,  (3.2) 

r  s  rV,.,  circulation 

w 

rQ  constant  in  characteristic  equation,  (3.2) 

At  t-^nterval  size  in  finite  difference  solution 

e  l/(k£Re) 

st  l/(kt  Re) 

t  =  (rV)p/(2r),  nondimensional  vorticlty 

G  azimuthal  angle 

<  constant  in  expression  for  radial  velocity  with  laminar 

Ekman  layer  (See  (2.8).) 

v  kinematic  viscosity  of  fluid 

p  density  of  fluid 

0  angular  velocity  of  spinning  cylinder 
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